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Abstract: We investigate the leptophilic properties of Dirac gauginos in an R- 
symmetric N=2 supersymmetric model with extended gauge and Higgs sectors. The 
annihilation of Dirac gauginos to leptons requires no chirality flip in the final states 
so that it is not suppressed as in the Majorana case. This implies that it can be 
sizable enough to explain the positron excess observed by the PAMELA experiment 
with moderate or no boost factors. When squark masses are heavy, the annihilation 
of Dirac gauginos to hadrons is controlled by their Higgsino fraction and is driven by 
the hZ and W + W~ final states. Moreover, at variance with the Majorana case, Dirac 
gauginos with a non- vanishing Higgsino fraction can also have a vector coupling with 
the Z gauge boson leading to a sizable spin-independent scattering cross section off 
nuclei. Saturating the current antiproton limit, we show that Dirac gauginos can 
leave a signal in direct detection experiments at the level of the sensitivity of dark 
matter searches at present and in the near future. 
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1. Introduction 



N=l supersymmetry broken around the weak scale would be the prime new physics 
candidate to be searched at the LHC. It provides not only an appealing explanation 
for the origin of the electroweak symmetry breaking but also a natural dark matter 
(DM) particle of the Universe. While N=l supersymmetry is enforced by the chiral 
structure of matter fields in the Standard Model, the gauge sector may be extended to 
have N=2 supersymmetry in which gauginos can be Dirac particles [JJ, 0, 0, |[ [|, ||. [7]. 
|8|. (for a review of N=2 supersymmetry, see |10j)- Determining the Majorana/Dirac 
nature of gauginos will be an interesting task for future experiments looking for 
supersymmetric CP and flavor violation |TT, [12|, collider signatures [ I3"| , [Tl| , and 
dark matter properties |0J [IB], |T7|, [T8| ]. 

There are some important distinctions between Majorana and Dirac gauginos. 
First, the annihilation of a Dirac gaugino pair into a fermion-antifermion pair, XX ~^ 
ff, has a non- vanishing s-wave contribution even in the limit of vanishing fermion 
masses, and thus the leptonic final states are not suppressed. This implies that a 
Dirac gaugino can provide a viable DM explanation for the abundance of energetic 



electrons and positrons recently observed in cosmic rays |19|, |5D, 7L, 22], ^3|, [24], . 
Second, Dirac gauginos can have a vector coupling with the Z gauge boson, xi^X ^ \ 
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leading to sizable spin-independent scattering off nuclei. As a consequence, the 



Higgsino component in the DM gets constrained by the direct detection data [26 



and also by the non-observation of antiproton excesses in PAMELA |?7| as will be 
discussed in detail in the following. 



In order to ensure such Dirac nature of DM, Majorana mass terms, which provide 
the mass splitting between the Dirac components, have to be highly suppressed. 
Otherwise, the heavier component of the two quasi-degenerate Majorana gauginos 
will decay to the lighter one, and the galactic DM will consist of a pure Majorana 
gaugino. If one assumes the N=2 structure in the Higgs sector, i.e. that the two 
Higgses H u and Ha form an N=2 hypermultiplet ]7|, the Dirac structure imposed 
at the tree level is spoiled by a large amount due to radiative corrections with Higgs- 
Higgsino and fermion-sfermion in the loop. A way to enforce an (almost) pure Dirac 
property for gauginos is then to assume a continuous R symmetry which forbids the 
usual /i term for the Higgs bilinear coupling H u Hd and the soft breaking trilinear A 
term [12fl. 



In this paper, we will work out a framework for R-symmetric Dirac gauginos, 
which are introduced by extending the gauge/Higgs sector of a supersymmetric 
model. We assume the Dirac bino as the main component of the lightest super- 
symmetric particle. The annihilation of two Dirac binos into a lepton-antilepton 
pair usually dominates over the annihilation into a quark-ant iquark pair, as slep- 
tons are typically much lighter than squarks. As we will see, the annihilation to a 
CP even final or intermediate state turns out to be velocity-suppressed, which is a 
property that is also true in the usual Majorana gaugino case. Then, the Dirac bino 
annihilation to the two W bosons and that to a CP even Higgs boson and 

a Z boson (hZ), both controlled by the Higgsino component in the DM particle, turn 
out to be the next important sources of the cosmic antiproton flux. Quantifying the 
suppression factor for the Higgsino component, one can make a prediction for the 
direct detection rate. 



In Section 2, an R-symmetric N=2 extension of the Minimal Supersymmetric 
Standard Model is introduced to enforce an almost pure Dirac gaugino structure. 
Based on this framework, we derive interaction vertices for the Dirac bino, which is 
assumed to be the dark matter, and then provide a qualitative analysis for its indirect 
and direct detection properties in Section 3. A numerical analysis is performed in 
Section 4 to obtain various fits and constraints of the leptophilic Dirac bino dark 
matter from the current indirect and direct detection data. We give our conclusions 
in Section 5. 
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2. The R— symmetric gauge/Higgs sector and Dirac bino cou- 
plings 



Let us first note that a small Majorana mass term, breaking the Dirac structure, will 
produce quasi-degenerate two Majorana gauginos, the heavier of which will then de- 
cay to the lighter one. Assuming a tiny splitting Stum in the Dirac bino components, 
one expects to have one-loop diagrams generating the magnetic moment operator: 



a' rrir, ~ 



4tt m? 



^B^ U B 2 F^, (2.1) 



between two almost degenerate Majorana components (m^ 2 — = Smu) where 
trig pa rrig i2 is the almost degenerate bino mass and m is a sfermion mass. The decay 
rate for the process B 2 — > B^ is 



a ( a' m| 
4tt \ 4tt m 2 



2 

5m M . (2.2) 



In order to maintain the Dirac gaugino structure for the present DM, we require that 
its lifetime is larger than the age of the Universe. This gives us 

5m M < l(T 33 GeV, (2.3) 

for nig ~ rh. Thus, the Majorana gaugino mass term has to be suppressed below the 
order of Stum "C m|/ 2 /M|>, where m 3 / 2 is the gravitino mass ~ 1 TeV. In N=l super- 
gravity, a nonzero value of the superpotential Wq is required to tune the cosmological 
constant. As wq breaks an R symmetry, it can usually generate Majorana gaugino 
masses of the order of m 3 / 2 ~ wq. In fact, avoiding such a leading contribution to 
gaugino masses and A or B soft supersymmetry breaking terms has been discussed 
in the context of split supersymmetry PBfl . Furthermore, a loop contribution from 



anomaly mediation was also shown to be further suppressed |2S[ Concerning the 
anomaly mediation contribution, it is interesting to remark that, saturating the lower 
bound, one finds 5rriM ~ (5 ,2 /167r 2 )(l/167r 2 )(m3 / , 2 /Mp) |28j which can be compatible 
with the condition (| 



An almost R-symmetric supersymmetry breaking can be realized in the following 



way. Let us introduce an extended gauge-Higgs sector [12 



W = fXiH d R u + fi 2 R d H u - V2g a (^H d T a R u + &H u T a R d ) $ a , (2.4) 

where $ a = (<f) a ,ifj a ) is the N=2 counterpart of the gauge superfield W£. In our 
phenomenological approach, the N=2 structure is supposed to be broken in an ap- 
propriate way leading to an effective N=l superpotential, as discussed below. For 
instance, £1^ in Eq. ( ^.4[ ) are arbitrary parameters breaking the N=2 relation. Note 
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that the usual Higgs superfields H^ u cannot form an N=2 hypermultiplet as the cor- 
responding fi and B terms, breaking an R symmetry, will generate a large Majorana 
mass mi/2 oc (a/47c)fiB at one-loop. A similar one-loop contribution arises also from 
the presence of left-right mixing sfermion masses, depending on the /i or A term. 
Thus, we are forced to extend the Higgs sector to avoid the usual /x and B terms 
and we will consider a situation in which the N=2 structure is not respected in the 
extended Higgs sector. Now let us introduce a supersymmetry breaking field X and 
two R symmetries, denoted by 1Z and J ', which are broken by two order parameters 
Xq and wq. The R charge assignment is given by 





W a a $ a H dtU R u4 


X 


Xq Wq 


K 


12 





2 2 


J 


10-13 


2 


2 



Note that the two R symmetries are not compatible with the N=2 structure [[Hj as 
different J charges are assigned to H d u and R u> d and thus they are not considered 
to come from the same hypermultiplets. The superpotential of the supersymmetry 
breaking sector can be written as 

W = w + x X } (2.6) 

which gives (X) = 9 2 Fx, where Fx = £o ~ ^3/2 ~ wo is assumed. Now one can see 
that the Dirac gaugino mass and the extended /1 terms can come from the following 
F or D terms: 

XW a a Q a $ a \ F1 HR$\ F or HRX\ , (2.7) 

where Q a is the supersymmetric covariant derivative. The last two terms in the 
above equation reproduce the superpotential ( |2.4p and the first term gives the Dirac 
gaugino mass term: C so f t = M a %p a W a . In the following, we will concentrate on the 
Dirac gauginos in the electroweak sector. The Dirac bino (wino) mass is denoted by 
Mi (M 2 ). The Majorana gaugino mass (and similarly the A term) gets a contribution 
from 

Wy/*WlXX^\ D . (2.8) 

This gives the (Majorana) soft term C so ft = 5rriMW a W a with bmu ~ m l/2/M P 
which can be compatible with the bound ( |2.3|) assuming a small coefficient of the 
order of 10~ 4 for m 3 / 2 ~ 300 GeV. Thus, we can safely neglect the Majorana mass 
terms in our discussion. 

Such a framework allows us to have highly suppressed Majorana mass terms and 
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results in the following R-symmetric gaugino-Higgsino mass matrix: 



M-- 






Mi 





-£im z s w c p 





£ 2 m z s w s /3 


Mi 





-m z s w Cf3 





m z s w Sf3 








-mzSwC/3 





/'i 








-iimzswcp 





/'i 














m z s w sp 













^mzswsp 











f>2 






(2.9) 



in the basis of , B, Hd, R u , H u , Rd). Note that, in order to make a transparent 
discussion of the Dirac gaugino properties, we reduced the number of free parameters 
by decoupling the wino component which can be justified in the limit of heavy wino 
mass: M 2 3> Mi, 

The Dirac gaugino mass matrix ( |2.9|) can be diagonalized by 6 angles which can 
be parameterized as 



ci si 

C 4 S 4 

-s 1 ci 00 

-s 4 c 4 

10 

1 



c 2 




s 2 





c 3 






s 2 
s 3 
10 
10 
c 2 



-s 3 c 3 



1 
1 
c 5 

-s 5 








c 6 








s 5 

s 6 

c 5 

-s 6 c 6 



(2-10) 



The six angles are defined by the diagonalization condition: Aidiag = N T M.N , 
where the diagonalized mass matrix Aidiag takes the diagonal form of 2 x 2 sub- 
matrices which have non-vanishing symmetric off-diagonal components of the Dirac 
mass eigenvalues denoted by m x , m^ Hi , and m^ H . 

The sizes of the mixing angles play an important role in determining the cosmic 
antiproton flux and the direct detection rate as we will discuss in detail in later 
sections. It is now instructive to find their approximate expression in the limit of 
small mixing angles: Sj -C 1. In the leading order of m z sw Mi < /ii i2 , the four 
mixing angles are given by 



«3 



S2 



cim z s w c p (M 1 + jinx) 
Ml - 

c 3 c 4 m z s w s^(c 2 1 ^ 2 M 1 + fj, 2 ) 



c 1 c 2 m z s w s p (clM 1 + £ 2 /x 2 ) 
c\c\Ml - & 



s 4 



and correspondingly the approximate mass eigenvalues are 



(2.11) 



m 



x ~ (ciC 2 Mi + s l c 2 m z s w cp - s 2 m z s w sp)c z c i 

- (cic 2 ^ 2 m z s w si3 - s 2 /i 2 )s 3 + (c 1 c 2 ^im z s w C/3 + sic 2 /ii)c 3 s 4 , 
m^ Hi « (siMi - cim z s H /c /3 )s 4 + (-s 1 £ 1 m z s w c f 3 + ci/ii)c 4 , 
m 4>H 2 ~ (S3C4M1 + c 3 ^ 2 m z s w s^)c 1 s 2 + (s 3 c 4 m z svyS/3 + c 3 /i 2 )c 2 . (2.12) 
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Also in the limit of small s« and = 1, the remaining two mixing angles have the 
form of 

„„ i ,„I!!l£s fa + lV (2.13) 

After the diagonalization, one can find the following interaction vertices of the 
dark matter particle X relevant for our analysis: 

C f = V2g'Y fR Af 22 [fP LX fR + xP R ff* R } 

+ V2g'Y fL M 22 [fP RX c f L + TPlHI] , (2.14) 

C h o = ^h°[C 1 ^ Hi P RX c +rPLi'H 1 ) 

+ C 2 $ H2 P rX c + TPl^h 2 ) 
+ C'S Hi P lX c + TPr^h 1 ) 

+ c> 2 (JP H2 Plx c + TPr4>h 2 ) 

+ 2S sXX ], (2.15) 
Cz = ^ Z ^H 1 l,PLX c + ri,P L ^H 1 ) 

+ Si{^h 1 1»PrX c + Ti^Pr^hJ 
+ ^h 2 1^PlX c + Ti^PlM 
+ ^2^h 2 1^PrX c + Ti^Pr^h 2 ) 
+ $vXl»X + ^aXI^x] , (2-16) 

where 



Ci = s a Af 22 J\f 33 + c a Af 22 Af 53 , C 2 = c a Af 22 Af 55 + s a J\f 22 J\f 35 , 

5s = ^ (sJVmA^i + C a Af 22 Af 51 + ilSaNuNv + 6CaMlAT 62 ) , 

5i = -M4M2 + A^A^ , 5 2 = +A" 66 A' 62 - A^A^ , 
5; = -A^A/^ + A'sgAk , 5 2 = +A/" 55 A/'5i - A^sA^i , 

^ = ^(^3 2 i-A^ 1+ A^ 2 -A4 2 2 ), 

5\ = i(A4 2 ! - A?, - A" 4 2 2 + A" 6 2 2 ) . (2.17) 

Here Yf L/R is the U(l) hypercharge of the left /right-handed fermion j^j R and the 
corresponding sfermion J'l/r- The angle a is the usual diagonalization angle of the 
neutral CP even Higgs bosons. In order to simplify the discussion in our analysis, 
we will decouple all heavy Higgs bosons by assuming ^> m z , where ttia is the 
pseudoscalar Higgs boson mass. 
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3. Indirect and direct signals of Dirac neutralino DM 



In the standard Majorana case, the direct annihilation of two neutralinos to two 
fermions requires a helicity flip in the final states and is vanishing in the limit of 
massless fermions, mf = 0. As a consequence of this, the annihilation cross section 
to light fermions is suppressed by a factor of m*/m^. On the other hand, when 
the neutralino is a Dirac particle this suppression is not present, so that the direct 
annihilation to leptons is largely enhanced and can be the dominant annihilation 
channel. In this case, the t-channel annihilation to quarks through squark exchange 



can be suppressed if m§ mj pO 



Then, in the presence of a large Higgsino component in the Dirac neutralino 
composition, annihilations to gauge and Higgs bosons can also lead to significant 
hadronic final states. Among them, annihilations to quarks through a Higgs particle 
exchange are suppressed either by a heavy pseudoscalar Higgs boson in the decoupling 
limit, or by a suppressed s-wave contribution in the non-relativistic limit, an effect 
which is present also in the Majorana case. 

Under these conditions, we analyze the dominant annihilation channels and the 
direct detection bound of the Dirac neutralino based on the results obtained in the 



previous section. In Sections |3.1| and [372], along with the dominant annihilation of the 
Dirac neutralino to leptons we will discuss the other relevant annihilation channels 
to hadronic final states. We will present the neutralino-nucleon cross sections in 
Section |3~3] . The result of this section will be used in Section [O] to constrain the 
Higgsino component of the Dirac neutralino. 

3.1 Annihilations to 11 and hZ 

First, let us consider the t-channel annihilation to leptons through slepton exchange 
which has to be the dominant annihilation channel to explain the observed energetic 
electron and positron spectrum. In the case of mr <C mj , we have 

<™>« - Y* , , ™S . (3.1) 



c 



w 



R (m~ + ml) 
Ir x 



On the other hand, the case of mr ~ mr — mr leads to 



Mflg<^*f ^ 2 m L^ (3-2) 



% 1 (mf + m2)2 

where = Yf + Yf . For the three families of the sleptons having same masses, the 
above expression must be multiplied by a factor of 3. 

The Dirac gaugino t-channel annihilation through Higgsino exchange produces 
a Higgs and a Z bosons, which can significantly contribute to the cosmic antiproton 
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flux. The annihilation rate for this mode is 



(av) 



hZ 



mi 



2 

mi 



5 r 4 2 



mi 



(3.3) 



2 mi 



i\2 



+ ^y + y') 2 



where 



E 



m x m^ H Ci5i 
mf 



E 



i 



2 

mf 



^ m 2 C^ 



2^ ' 



with m\ = + m 2 x — m|, m| = (m 2 z + m^)/2. Here — m| <C m 2 is assumed. 

In Section |4], we will present numerical results on the limit of (av)hz coming 
from cosmic antiproton fluxes. It is, however, useful to see a qualitative behavior of 
the small mixing limit: 5i, 8[ <C 1. Taking = 1 in Eq. ( p. 11 ), the cross section (|3T 
reduces to a simple approximated form in the zeroth order of m\: 



(av) 



mr 



hZ 



ir a 



4c^ (m| + m 2 ) 2 



m 



H 



sin 2 (a + ft) 



3D becomes 



where m^ H ^ ~ ^t/> H2 — wig and Mi <C is assumed. The ratio between 

1 (m| + m 2 ) 2 



hZ 



81? 



(<™}fl 

Assuming now mr ~ mg and Y z 4 

(crv)hz 



m~ + m 2 ) 2 



mg 



sin (a + ft) . 



(3.4) 



and 



(3.5) 



51/16 (three even leptophilic), we have: 



< 



51 



i _ UhL 



m H 



< 



51 



(3.6) 



Comparing the annihilation cross section to charged leptons compatible with the 



PAMELA positron excess [21] and the upper bound on the annihilation cross section 
to hZ from the PAMELA antiproton data |27| , presented in Fig. ^| of Section |], we 
obtain a conservative limit of (av )hz/(°~v) l - l < 0.3 which can be easily satisfied, as 
shown in Eq. ( |3.6| ). In addition, even for only one lepton channel, the condition is 
also satisfied: (av)hz/ (& v )u < 2/17 < 0.3. 

3.2 Annihilations through Z 

In the zero-velocity limit, the following annihilation channels through s-channel 
exchange of a Z boson have non- vanishing amplitudes. Compared to Eq. ( |3.3| ), these 
rates contain an extra factor ~ (^\l^\w) $ 2 anc ^ thus can be important in certain 
parameter regions. 
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The annihilation rate for XX ~^ ff is 



(av) 



ff 



2na 2 

,4 „4 

w 

2 



m: 



m, 



c w s w 



(3.7) 



m. 



rrir 



-qWa + 



m: 



(Am 2 



m 



2 \2 



2(g 2 v + g 2 A ) + ^(g 2 v -2g 2 A ) 
m x 



£ 4 



where gy = (T 3 — 2QfS 2 v )/2 and = — T 3 /2. The annihilation rates for — > hZ 
and W + W~ are respectively given by 



(<7U) 



ftZ 



(av) 



WW 



■na 2 \ 


/ m| — m| 


4 4 






m x 




-<) A ~ 






na 2 m x 


\j m l- m w 


s w ( 4m x ~ m D 2 



(3. 



m; 



(Am 2 , 



w 



2 \2 



1 + 



2m? 



m 



m; 
+ 16^ 



mr 



rf 4 



17-3- 



m 



w 



rn 



m 



(3.9) 



Finally, there is also an annihilation channel XX ~^ 2/i. The annihilation rate for this 
channel is given by 



(av) 



na 



mr 



m; 



;) 3/2 



2/( 



4r 4 <? 4 



2 °V 



(3.10) 



"W°w m x (Am 2 x 

Among these four channels, the process XX ~~ >* %h is the smallest one. The hZ channel 
is larger than the // channel even for the top quark. Thus, the dominant channels are 
hZ or W + W~. The ratio between the two channels is approximately 5 A : Sy/A. The 
ratio between the cross sections for annihilations to 11 and W + W~ (hZ) depends 
on two parameters, m x and 5y ($a)- The upper bound on the ratio between the 
annihilation rates for XX ~~ M an d W + W~ (hZ) is approximately 0.5 (0.3), as can 
be seen from Fig. ^|. Then, comparing Q3.9|) and (|3.8|) with ( |3.2D , one finds the limits: 
<V ^ 0.3 — 0.03 and 5 A < 0.2 — 0.02 depending on m x in the mass range 200 GeV 
< m x < 2000 GeV. Note that the upper limit of 5y is close to the current sensitivity 
of the direct detection of DM as will be discussed in the following subsection (see 
Fig. [[]). This will be also confirmed by our numerical analysis performed in a more 
general parameter space. 

3.3 Direct detection through Z and h exchange 

The vector interaction via t-channel Z boson exchange (|2.16 ) leads to a spin- 
independent neutralino-nucleon cross section. The cross section for the neutralino- 
nucleon vector interaction is given by 

1 -l 2 



a. 



vector 



4 



na 



64s 



w c w 



mi 



2s 



w) 



(3.11) 



- 9 - 



' CDMS II 


<! v =0.3 — ' 






(S s =0.8 


1.2 


S s =03 


(! v =0.1 






Ss =0.1 


(S v =0.07 



200 400 600 800 1000 

(GeV) 

Figure 1: Exclusion plot for the spin-independent neutralino-nucleon cross section asi- 
The (blue) thick solid lines are the cross sections via the vector interaction corresponding 
to 8y = 0.3, 0.2, 0.1, and 0.07, respectively. The (red) dashed lines are the cross sections 
through the Higgs exchange corresponding to Ss = 0.8, 0.3, and 0.1, respectively, taking 
m h = 115 GeV. The dot-dashed line shows the CDMS II limit @. 



where /i n is the reduced mass for x-nucleon, and Z and A are the atomic number 
and the atomic weight of the target nucleus, respectively. 

In addition to the vector interaction, the Dirac gaugino x interacts with nucleons 
through t-channel Higgs exchange fl2.15|) . The neutralino-nucleon cross section for 



the scalar interaction is 



p 2 0.14 2 x Ag' 2 h 2 hss filmj 

^scalar ~ °S ~ ^4^2 ' V 3 - 1 ^) 



7T m^mi 



where hhss is the Higgs-s quark-s quark Yukawa coupling. In the decoupling limit, 
this result is estimated as 



0.14 2 x 16vra 2 fi 



2 ml 



^scalar - °S 2 _2 ^4^2 " 

b W c W " L h" L W 

In Fig. [I], we present the neutralino-nucleon cross sections via the vector inter- 
action (thick solid lines) as a function of m x for Sv = 0.3, 0.2, 0.1, and 0.07. The 
scalar interaction cross sections (dashed lines) are also shown as a function of m x for 
5s = 0.8, 0.3, and 0.1. The limit from the CDMS II experiment is shown as a 
dot-dashed line. 
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4. Analysis 



In this section, we analyze the phenomenology of Dirac gauginos in a quantitative 
way. As discussed in the previous sections, Dirac gauginos can easily meet the 
requirements of a leptophilic DM candidate, and so we will start in Section |4.1| by 
summarizing present constraints on generic leptophilic models. We will then proceed 



in Section [4.2| to discuss specific properties of Dirac binos as leptophilic DM, while 



in Section JO we will extend the analysis to Dirac neutralinos of generic composition 



in some particular examples. 

4.1 Constraints on leptophilic dark matter 

Leptophilic DM annihilation is usually advocated in order to explain simultaneously 
the PAMELA positron excess |H] and the excellent agreement between the observed 



antiproton spectrum and the corresponding standard expectation [p7| . The present 
situation of a leptophilic DM candidate annihilating democratically into charged 
leptons of the three families is summarized in Fig. [|, where the annihilation cross 
section (erf) is plotted as a function of the DM mass m x . In particular, the thick and 
thin closed solid contours show the range of values compatible with the PAMELA 
positron excess [[Zl|] and the FERMI-LAT e + + e" data |[25|| , respectively. The con- 



tours are calculated by requiring that the \ 2 P er degree of freedom for fits to the 
data is less than or equal to 1 (only PAMELA data above 10 GeV have been included 
in the fit), and assuming for the secondary positron and electron backgrounds the 
conventional GALPROP model denominated by Model #0 in Table 1 of |JT|. As 



far as the propagation model for primary electrons and positrons is concerned, we 
have adopted a Navarro-Frenk-White (NFW) profile [^j with medium ranges for 
the diffusion coefficient and for the size of the propagation region, corresponding to 
the model NFW-med in Table 2 of [j33|| . The two dot-dashed lines in the upper part 



of Fig. represent conservative 2 a C.L. upper bounds for (cry) obtained from the 
flux of e + + e~ observed by FERMI [[Jl| (thin line) and the e + + e~ flux measured 
by HESS |24| (thick line). 

In the same Figure, we also plot for reference with the short dash-long dash 
line the upper bound on (av) obtained by considering the imprint on the Cosmic 
Microwave Background Radiation (CMB) from the injection of charged leptons from 
DM annihilations at the recombination epoch fl3"I , |55| , [J?J. In particular, the plotted 



line is obtained by taking the WMAP5 constraint on the quantity f(crv) from Figure 
4 of Ref. p4| , where / is defined as the average fraction of the DM rest energy 



deposited in the gas at z ~ 800-1000, and taking / ~ 0.39 from the analysis of 
Ref. p5| . As one can see, at present the effect of the CMB constraint on leptophilic 



models which can explain the PAMELA excess is limited to very large masses, m x > 
1.5 TeV. Future polarization data from Plank will be able to improve this kind of 
limits considerably [3^ . 
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Figure 2: Annihilation cross section times velocity {crv) for a DM candidate annihilating 
democratically to charged leptons of the three families as a function of the mass m x . 
The thick solid contour shows the range of values compatible with the PAMELA positron 
excess [21|; the thin solid contour is the range compatible to the observed e + + e~ flux 



measured by FERMI-LAT [25]; the thin dot-dashed line is the 2 a upper bound from the 
observed e + + e~ flux of FERMI; the thick dot-dashed line is the 2 a upper bound from 
the e + + e~ HESS measurement [24]; the short dash-long dash line represents the upper 



bound on {crv) from CMB [34, [35|, 36]. The lower and upper dotted lines show the upper 



bound on (cry) obtained by comparing the expected gamma-ray flux produced by Inverse 
Compton (IC) scattering of the final state leptons to the FERMI-LAT measurement of 
the diffuse gamma ray emission with or without subtraction of the expected standard 
background, respectively [p7| . The short dashed and long dashed lines respectively show 
the upper bounds on the annihilation cross section to the final states W + W~ and hZ from 
the PAMELA antiproton data |27|]. See Section | for details. 



Finally 1 , the dotted lines show the upper bounds on {crv) obtained by compar- 
ing the expected gamma-ray flux produced by Inverse Compton (IC) scattering of 
the final state leptons to the diffuse flux of gamma-rays measured by FERMI at 



1 Signals from DM annihilation in the Galactic Center are very sensitive to the choice of density 
profile. Therefore, we do not consider them here. 
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intermediate Galactic latitudes p7| . We do not include in our analysis the prelim- 
inary data points very recently released by the FERMI-LAT Collaboration. Due 
to the large experimental errors, they would not change our conclusions ||38|| . The 
upper dotted line has been obtained by directly comparing the flux expected from 
leptophilic DM (calculated in the approximation of Ref. |]39]]) to the FERMI-LAT 
data, while the lower dotted curve has been obtained by subtracting the estimation 
of the standard background from the data (as a conservative estimation of the latter 
we have assumed the lower boundary of the dashed region indicated as "Total" in 
the left-hand panel of Fig. 1 of Ref. [p7|| ). 

The IC analysis is described in detail in Fig. [| where the diffuse FERMI-LAT 
data at latitudes 10° < b < 20° from [37| as a function of the gamma-ray energy E 1 



are compared with the corresponding fluxes from leptophilic DM annihilations. In 
particular, the lower solid curves represent the gamma-ray fluxes from DM annihi- 
lation for the representative values m x =200, 500, 1000, 1500, and 2000 GeV (from 
bottom to top), where for each m x the value of (av) is taken from the lower bound- 
ary of the thick solid contour of Fig. [2] (that we adopt here and in the following as 
a conservative estimation for the annihilation cross section required to explain the 
PAMELA data). The same expected fluxes are also summed to the standard back- 
ground, shown by the thick solid line in the range 0.1 GeV < E 1 < 10 GeV, and 
the results are plotted in the same energy range. The dashed curves correspond to 
the predictions for the additional value m x ~ 1135 which, as shown in the magnified 
inset at the upper-right, can be adopted as an estimation of the upper bound for m x 
from the IC flux. 

A comment is in order here. Following other similar analysis in the literature 
[y], |33| , we take the results shown in Fig. || as only indicative of the present trends. 



In fact, it is clear that a consistent \ 2 analysis would require a robust assessment of 
the systematic uncertainties, an established knowledge of the expected backgrounds 2 
and a full marginalization over the astrophysical parameters, something which is not 
possible in light of the present uncertainties. 

4.2 The Dirac bino as a leptophilic DM candidate 

Dirac gauginos can easily meet the requirements of a leptophilic DM candidate, 
since their annihilation cross section to leptons can be dominant, not being chirality 
flip suppressed as in the case of Majorana gauginos, while annihilations to quarks 
potentially capable of producing an excess in the antiproton signal can be inhibited by 
assuming large masses for the squarks |3"0 |. If the additional hadronic annihilation 



channels discussed in Sections |0| and |37^ are further suppressed by assuming a 
negligible Higgsino component in the neutralino composition (by taking 3> M\ 



2 At variance with antiprotons, the expected standard background for electrons and positrons 
could in principle be only in part of secondary origin, with a sizable primary fraction (see for 



instance g(| 
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Figure 3: Diffuse gamma-ray flux produced by Inverse Compton scattering of the final 
state electrons produced in DM annihilations as a function of the gamma-ray energy 
for a leptophilic DM candidate annihilating democratically to charged leptons of the three 
families. The lower solid curves represent the gamma-ray fluxes from DM annihilation for 
the representative values m x =200, 500, 1000, 1500, and 2000 GeV from bottom to top, 
where for each m x the value of (av) is taken from the lower boundary of the thick solid 
contour of Fig. ^| (adopted as a conservative estimation for the annihilation cross section 
required to explain the PAMELA data). The same expected fluxes are also summed to the 
standard background, shown by the thick solid line in the range 0.1 GeV < EL < 10 GeV, 
and the results are plotted in the same energy range. The dashed curves correspond to the 
IC predictions for the additional value m x ~ 1135 which, as shown in the magnified inset 
at the upper-right, can be adopted as an estimation of the upper bound for m x from the 
IC flux from our analysis. The experimental points are the FERMI-LAT data at latitudes 



10° < b < 20° from J37| 



in the mass matrix of Eq. Q2.9|) ), a totally leptophilic model is obtained for which 
only the constraint m x < 1 TeV from the IC flux discussed in Fig. |3] is applied 3 . 
The phenomenology of such a DM candidate is quite simple, and is summarized in 



3 Notice that the IC gamma-ray flux retains the directional information of the source, allowing in 
principle to establish whether the e ± PAMELA excess is diffused in all the Galactic halo or localized, 
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Figure 4: Boost factors for a Dirac gaugino DM candidate annihilating democratically 
to leptons of the three families as a function of the mass and in the limit \fx\ ^> M\. 
Solid lines show the ratio between the annihilation cross section required to explain the 
PAMELA excess and the value of (av) calculated using Eq. ( |3.2| ), Dashed lines show the 
ratio between the observed lower value of the DM relic density $l m in 

h 2 = 0.098 H and the 

corresponding thermal relic abundance for the Dirac gaugino. Thin lines refer to mj/m x = 
1 (where we have assumed a small mass splitting between mj and m x sufficient to avoid 



coannihilations between sleptons and neutralinos) while thick ones to mj/m x 
mj is the slepton mass. 



-2, where 



Fig. where the solid lines show the ratio between the annihilation cross section (av) 
required to explain the PAMELA data (taken as before from the lower boundary of 
the thick closed solid contour in Fig. ||) and that calculated by making use of Eq. (|3.2| ) 
for mj/m x =1 (thin line) and 2 (thick one). In particular, mj/m x >1 maximizes the 
expected annihilation cross section by assuming the lightest possible slepton mass 
compatible with the requirement that the neutralino is the Lightest Supersymmetric 
Particle (LSP). In this latter case, the Dirac gaugino annihilation cross section to 
leptons can explain the PAMELA data for m x ~ 330 GeV, while m x ~ 200 GeV is 
needed if mj/m x =2. 

It is well known, however, that the values of the annihilation cross section that 

for instance, in a DM clump close to the Solar System. The preliminary data of Ref. |}7j do not 
provide this piece of information since they are averaged over all galactic longitudes 0° < I < 180°. 
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can explain the PAMELA excess ((crv) ~ 10~ 24 cm 3 s _1 ) are about two orders of 
magnitude larger than the range compatible to a standard thermal relic abundance in 
agreement with observation (2-3xl0~ 26 cm 3 s _1 at the decoupling temperature T ~ 
m x /20). This is naively shown in Fig. ^, where the dashed lines represent the ratio 
between the 2 a lower value of the observed DM relic density Q m i n h 2 = 0.098 p2[ and 



the corresponding thermal relic abundance for the Dirac gaugino. The thin dashed 
line refers to mj/m x = 1 while the thick one to mj/m x = 2. From this figure, one 
can see that a light Dirac gaugino (m x < 500 GeV) can explain the PAMELA data 
with moderate or no boost factor at all if, for instance, some non-standard expansion 



history of the Universe |43[ or non-thermal production mechanism is advocated in 
order to reconcile its DM relic abundance to the observation. Alternatively, a thermal 
relic density in agreement with observation is attained for m x > 800 GeV, but in this 
case a large enhancement of the annihilation cross section to leptons is necessary. 
Such an enhancement could be provided by an astrophysical boost factor due to 
substructures, although recent analysis seems to disfavor values larger than order 10 



HJ]. Furthermore, this piece of information can be combined to the data on e + + e~ 
from FERMI-LAT, which, as shown by the thin closed contour plot in Fig. ||], imply 
a lower bound on the neutralino mass m v > 465 GeV. Notice, however, that lower 
values for the mass of the annihilating DM particle can be in principle assumed 
if some additional contribution to the electron+positron background is claimed at 
energies E > 1 TeV to explain the FERMI-LAT data independently j40j p . 



4.3 The case of a Dirac neutralino of general composition 

As discussed in Section p.3|, an important phenomenological feature of a Dirac neu- 
tralino is the presence of a vector coupling with the Z boson, which on the other 
hand is vanishing in the Majorana case. This implies that, while in the standard 
Majorana scenario a diagram with the Z boson exchange can only contribute to the 
spin-dependent neutralino-nucleon cross section, 4 for the Dirac case it can lead to 
a much more sizeable cross section enhanced by a factor of ~ A 2 (see Eq. ( p.lip ). 
However, note that the coupling between a neutralino and a Z boson vanishes (as 
long as the coupling to the Higgs bosons) in the limit |//i,2| S> Mi of a pure gaugino. 
This implies that the detection of Dirac neutralinos through direct searches is possi- 
ble only in the presence of a Higgsino component in a Dirac neutralino of arbitrary 
composition. In order to explore this possibility, we discuss in this section the full 
parameter space of our model, by allowing | /Xx,2 1 — M\, and discuss the constraints 
coming from the hadronic annihilation channels discussed in Sections |3J] and |3.2| . In 
particular, in Figs. |5] and ^|, the M\ — fii parameter space of the Dirac neutralino is 
explored for the two representative cases /ii//i2 = 1 and = 0.1. Our numerical 

analysis is performed in the case of £i = £2 = 1. In each figure, the left-hand 

4 The most stringent limits on pure WIMP-proton and WIMP-neutron spin-dependent cross 
sections are given by KIMS El and ZEPLIN-III El , respectively. 
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Figure 5: Parameter space for a Dirac neutralino of general composition, for Hi=H2 
and tan/3=3,10 (left and right panels, respectively). In each panel, the region below the 
slanted solid line is excluded by the PAMELA antiproton data, while along the two dotted 
lines the ratio between the neutralino-nucleon scattering cross section calculated with Eqs. 
[( 3.11[ ), Q3.12D ] and the experimental limit from CDMS II is fixed to 1 (lower curve) 
and 10 -3 (upper curve). The regions below the lower dotted curves, where this ratio is 
larger than 1, are excluded by the direct detection constraint. Moreover, in each plane, 
the vertical thin dashed line to the left represents points where the boost factor, defined 
as the ratio between the annihilation cross section (cry) required to explain the PAMELA 



excess and that calculated by making use of Eq. (3.2), is equal to 10, while the vertical 



thick dot-dashed line to the right represents points at constant mass m x =l TeV. 



panel corresponds to tan/3 = 3 and the right-hand panel to tan/3 = 10. In all the 
Mi— fa planes, the predominantly Higgsino region below the slanted solid line is ex- 
cluded by an antiproton flux above the PAMELA data of Ref. [p7|, driven by the 
cross section to W + W~ of Eq. Q3.9Q and a subdominant contribution from the hZ 
final state of Eq. ( |3.3j ) (in the calculation of the latter cross section we have assumed 
a super symmetric Higgs sector in the decoupling limit, inn, a 3> m^, with m/i=115 
GeV). The corresponding upper bounds to these two cross sections, which as already 
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Figure 6: The same as in Fig. |5], with /X2//xi=10. In both panels, the three dotted 
lines show points where the ratio between the neutralino-nucleon scattering cross section 
calculated with Eqs. [( p. 11 ), ( |3.12| )] and the experimental limit from CDMS II |26|| is fixed 
to 1 (lower curve, barely visible in the right-hand panel), 1CP 3 (middle curve), and 10 -5 
(upper curve). 



mentioned in Sections |3J] and |3.2| turn out to be the dominant hadronic ones in 
our R-symmetric scenario, are plotted for reference in Fig. || as short-dashed and 
long-dashed lines, respectively. They have been obtained by conservatively requir- 
ing x 2 < 40 when comparing the theoretical expectation to the PAMELA data (this 
value corresponds approximately to a 99.5 % C.L. for 17 degrees of freedom). For 
the calculation of the antiproton fluxes, we have adopted the NFW-med model for 



propagation |33| and the parametrization of Ref. |48[ for the secondary antiproton 
background. 

In Figs. [5] and |6], the vertical thin dashed line to the left represents points where 
the boost factor, defined as the ratio between the annihilation cross section (av) re- 
quired to explain the PAMELA excess and that calculated by making use of Eq. ( |3.2|) , 
is equal to 10. In particular, the points to the right of this vertical line require a larger 
value. Note that an astrophysical enhancement to the effective cross section provided 
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by DM substructure larger than order 10 appears to be unlikely Moreover, in 
each M\-\ii plane, the vertical thick dot-dashed line to the right represents points 
with a constant Dirac neutralino mass rn x —l TeV. Configurations to the right of this 
line have larger masses and appear to be disfavored by the IC flux shown in Fig. [3[ 
Note that this is the only irreducible astrophysical bound which is constraining also 
in the pure gaugino limit 3> Mi. 

In the regions of Figs. |5] and ^| that are excluded by antiproton fluxes (below 
the slanted solid lines), a sizeable Higgsino fraction can also drive the neutralino- 
nucleon cross section above the present experimental limit. This is shown by the 
dotted lines which represent configurations with a fixed ratio between the expected 
neutralino-nucleus cross section and the present experimental upper bound from 
CDMS II p6|. In particular, in all the /ij — M\ planes this ratio is equal to 1, 10~ 3 , 
and 10~ 5 (if present), starting from the lower curve. In the regions below the lower 
dotted curves, this ratio is larger than one, so that they are excluded by the present 
limit from CDMS II. 

The expected fluxes for positrons, electrons, and antiprotons, as well as the 
expected direct detection cross section for a Dirac neutralino are shown in detail in 
Fig. |7] for the specific example m x ~ 465 and (av) ~ 1.4 x 10 -24 cm -3 s _1 . This 
configuration corresponds to the lower mass bound implied by the FERMI-LAT 
e + + e~ data if no additional sources of electrons and positrons are assumed and, as 
shown in Fig. || implies a boost factor of around 3 in the case of mj/m x ~1 and in 
the limit of a pure bino. In particular, in Figs, ^a and 0b, the corresponding positron 
and e + + e~ fluxes are plotted as a function of the energy, while the antiproton flux 
and the neutralino-nucleon cross section are plotted in Figs. 0c and [TH, respectively, 
for the case of a Dirac neutralino of general gaugino-Higgsino composition in the 
representative case of tan (3 = 3 and /ii//i2=l. Note that in the case of a non- 
negligible Higgsino-gaugino mixing the boost factor for the leptonic annihilation 
cross section required to explain the PAMELA data, which is plotted in Fig. f| in 
the pure gaugino limit 3> Mi, can be somewhat larger due to the suppression 
of the bino coefficient A/22- However, this effect is not sizeable. For instance, for the 
maximal Higgsino composition compatible to the CDMS II limit shown by a star in 
Fig. |7|d the boost factor turns out to be ~4 instead of ~3. 



5. Conclusions 

Determining the Majorana/Dirac nature of gauginos will be an interesting task for 
future experiments to look for supersymmetry. Moreover, as an interesting variance 
of the standard Majorana case, Dirac gauginos can be natural realizations of lep- 
tophilic DM, which has been recently proposed to explain the rising positron flux 
observed in cosmic rays by PAMELA without producing excesses in the antiproton 
signal. In this paper, we have analyzed the phenomenology of Dirac gauginos in 
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Figure 7: Various observables calculated for the representative configuration m x = 465 
GeV and (cu),t = 1.4 x 10 -24 cm 3 s _1 . (a) Solid line: positron fraction compared to the 
PAMELA data |l|. Dashes: secondary background (Model #0 of Ref. @); (b) Solid 
line: electron+positron flux compared to the FERMI-LAT data Dashes: secondary 
background (Model #0 of Ref. [|3l|j ); (c) Solid line: antiproton fraction compared to the 
PAMELA data p7] . Dashes: secondary background |4lj|] . The value of the hadronic anni- 
hilation cross section (av) w + w - used for the calculation corresponds to the configuration 
indicated with a star in panel (d) which is at the level of the experimental CDMS II exclu- 
sion plot; (d) Neutralino-nucleon cross section as a function of [i\ = \X2 for tan (3 = 3. The 
horizontal line is the corresponding CDMS II upper bound [Eg]. 



a specific supersymmetric realization, where a continuous R symmetry is assumed 
to protect vanishingly small Majorana masses. We have shown that a light Dirac 
gaugino (m x < 500 GeV) can explain the PAMELA data with moderate or no boost 
factor for light (> m x ) slepton masses. Furthermore, in the case of a non-vanishing 
Higgsino fraction, Dirac gauginos can have a vector coupling with the Z gauge boson 
leading to a sizable spin-independent scattering off nuclei. In some specific examples, 
we have shown that present constraints on the Higgsino fraction of Dirac gauginos 
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from direct detection experiments are at the same level of those coming from an- 
tiproton data. This implies that a Dirac gaugino signal is potentially at the level of 
the sensitivity of direct detection experiments at present and in the near future. 
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